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Nonlinear energy transfer and dissipation in Alfven wave turbulence are analyzed in the first 
gyrokinetic simulation spanning all scales from the tail of the MHD range to the electron gyroradius 
scale. For typical solar wind parameters at 1AU, about 30% of the nonlinear energy transfer close 
to the electron gyroradius scale is mediated by modes in the tail of the MHD cascade. Collisional 
dissipation occurs across the entire kinetic range k^pi > 1. Both mechanisms thus act on multiple 
coupled scales, which have to be retained for a comprehensive picture of the dissipation range in 
Alfvenic turbulence. 


Introduction. Spacecraft measurements find a radial 
temperature profile of the solar wind which can only be 
explained by the presence of heating throughout the he¬ 
liosphere [T]. The key mechanism of heating in the in¬ 
ner heliosphere up to ^20 AU is thought to be the dis¬ 
sipation of turbulent fluctuation energy, and its under¬ 
standing and description is one of the outstanding open 
issues in space physics [2]. Over the past decade, nu¬ 
merous studies, both observational |3H2] and theoreti¬ 
cal/ computational EHni, have focused on this topic, ex¬ 
tracting ever more sophisticated measurements of solar 
wind fluctuation properties, and accomplishing increas¬ 
ingly detailed turbulence simulations. 

As the solar wind plasma is only weakly collisional, 
a variety of kinetic effects such as cyclotron damping, 
Landau and transit time damping, finite Larmor radius 
effects, stochastic heating, or particle acceleration at re¬ 
connection sites can contribute to the conversion of field 
energy to particle energy, and thus determine how col¬ 
lisional dissipation will ultimately set in. A kinetic de¬ 
scription is crucial in order to judge the relative impor¬ 
tance of each of those effects. Due to the complexity of 
a nonlinear kinetic system, numerical simulations are es¬ 
sential to interpret observations and provide guidance for 
analytical theory. 

In the present Letter, we employ an approach based 
on gyrokinetic (GK) theory [TS], which is a rigorous 
limit of kinetic theory in strongly magnetized plasmas. 
Due to the assumptions of low frequencies (compared to 
the ion cyclotron frequency) and small fluctuation levels, 
the gyrokinetic model excludes cyclotron resonances and 
stochastic heating. In absence of these effects, we focus 
on the energetic properties of kinetic Alfven wave (KAW) 


turbulence, which has been demonstrated to be a crucial 
ingredient of solar wind turbulence mi- 

We address the following key questions: (1) Which 
spectral features can be found in a comprehensive simu¬ 
lation extending from the magnetohydrodynamic (MHD) 
range down to the electron gyroradius scale? (2) What 
are the characteristics of nonlinear energy transfer from 
large to small scales? (3) How is energy dissipated, and 
how is the dissipated energy partitioned between ions and 
electrons? 

Simulation setup. The nonlinear GK system of equa¬ 
tions is solved using the Eulerian code GENE [T7] to 
study the dynamics of KAW turbulence in three spa¬ 
tial dimensions. In order to model the energy injection 
at the outer scales of the system, a magnetic antenna 
potential, whose amplitude is evolved in time according 
to a Langevin equation |18| . is externally prescribed at 
the largest scales of the simulation domain. The driven 
modes are (0,1, ±1) and (1,0, ±1), where {i,j,k) are 
multiples of the lowest wave numbers in {kx, ky, kz), re¬ 
spectively. The mean antenna frequency is chosen to be 
Wq = O.Qujao (wao being the frequency of the slowest 
Alfven wave in the system), the decorrelation rate is set 
to 7 o = 0.7a;AO) and the normalized antenna amplitude 
is set to A|| 0 = coaoBoVS/C 2 k^Q'/NvA (setting 6 = 2, 
W = 4, C 2 = 1), in accordance with the critical balance 
condition a;iin a;ni [E]. 

The physical parameters are chosen to be similar to 
solar wind conditions at 1 AU, with /3i = SimiTi/= 1, 
Ti/Tf, = 1. Proton and electron species are included 
with their real mass ratio of mi!me = 1836. The elec¬ 
tron collisionality is chosen to be Ve = O.OOojao (with 
G = '\An/7m/r'e), a value small enough to not inhibit 
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kinetic effects, but large enough to reduce resolution re¬ 
quirements in velocity space. 

In order to maximize the effective dynamic range, the 
simulation domain is extended significantly compared to 
previously published work, to include scales larger than 
the ion gyroradius, allowing for a free distribution of en¬ 
ergy into the KAW or the ion entropy cascade m as the 
ion gyroradius scale is passed. The evolution of the gyro- 
center distribution is tracked on a grid with the resolution 
{n^,ny,n^,ny^^,n^,na) = (512,512,96,48,15,2). The 
plane perpendicular to the background magnetic field is 
resolved by 512^ fully dealiased grid points, covering a 
perpendicular wavenumber range 0.2 < kj_pi < 51.2 (or 
0.0047 < k±pe < 1-19), thus extending into the regime 
where electron finite-Larmor-radius effects become im¬ 
portant. Here, p^ = eB with the species index 

cr. The number of grid points in the perpendicular plane 
is thus increased by a factor of 36 with respect to the 
largest runs of this kind published to date |9]. 96 points 
are used to resolve the dynamics along the background 
field (the z direction), and 48 x 15 gridpoints are chosen 
to represent the (u||, /i) domain, where uy is the velocity 
along the guide field, and p = mv\/2BQ is the magnetic 
moment with respect to the guide field. The domain sizes 
in velocity space are chosen to extend up to 3 thermal ve¬ 
locities VTcr in both parallel and perpendicular velocities 
for each species a, where VTa = \/ 2T'o-/Wo-. 

Our simulations are performed using the same iterative 
expansion scheme as in Ref. [1], where simulations are 
initially run with low resolution and are then restarted 
several times with an increasingly fine grid, until the tar¬ 
get resolution is reached. The total runtime is chosen to 
span several antenna oscillation periods ta (in this case 
fend = 7.20 ta) in order to ensure that a quasi-steady 
state has been reached. 

Diagnostic methods. The key results of this study are 
obtained using a set of sophisticated energy diagnostics 
(partially introduced in Refs. [20H23]), which enable stud¬ 
ies of energy source, transfer and dissipation spectra sep¬ 
arately for each species, and which are applied to KAW 
turbulence for the first time here. In particular, we an¬ 
alyze the time derivative of the spatially averaged free 
energy density, which can be expressed in the case of an 
antenna-driven electromagnetic system as 

+ ~^C^||ant,(Tfc^ dtPak^ 

+ • (1) 

Here, the sum over k denotes a summation over all 
wavenumber pairs {kx,ky), and the angle brackets in¬ 
dicate a spatial average along the guide field. 
is the perturbed gyrocenter distribution, and hak = 


fak + {qa4>iak + tJ-Bi\\crk) is its nonadiabatic 

part. The overbar denotes an average over the gyro-ring, 
and Fqu is a Maxwellian background distribution with 
background density and temperature Tq„. The mag¬ 
netic potential Aiytoj ^ + Ayant.fc is understood 

to contain also the contribution due to the Langevin an¬ 
tenna Ayant.fc) which is necessary for a complete account 
of the energy contained in the system. The time deriva¬ 
tive dtPak = dt{fak + qaV\\Ai\\kFo„lcTQ„) is the quantity 
explicitly evolved in the GK Vlasov equation as imple¬ 
mented in GENE, and 


C = k^ 


STT^qjBo 


^ll'^o ('^fr) Eofjd'Uy d/r 




is a factor arising from the antenna-modified Ampere’s 
law, with Ao- = kj_yj2map/B^ql. By replacing dtPak in 
Eq. Q with any of the various terms contributing to its 
evolution, we can assess the impact of that term on the 
evolution of the free energy density. The nonlinear trans¬ 
fer function (i.e. the free energy balance contribution of 
the nonlinear term) thus reads 


Fkpq — 
X 




Tqo- 

Foa 


dv\\Ap[pa:qy - PyQa] [XlapKq 
+ qaV\\CA\\ 


Xlaq^<yp\ 


( 2 ) 


with k -\- p q = 0. Gompared to the definition used in 
Refs. [22l [2^. there is an additional term involving the 
antenna potential, and the electrostatic approximation 
has been dropped by using the full electromagnetic po¬ 
tential Xia = - '(^ll^illtot.^/c + pBi\\^/qcr. Note that 

the new antenna potential term does not satisfy the same 
symmetry properties as the rest of the transfer function, 
consistent with the fact that the antenna acts as an en¬ 
ergy source through the nonlinear term (but also through 
the parallel advection term). This source can be quanti¬ 
fied by measuring the symmetric part of the above trans¬ 
fer function. 

Field energy spectra. Before focusing on the nonlinear 
transfer physics, we analyze the spectra of the magnetic 
and electric field energy, which can be directly compared 
to spacecraft observations. As is common practice, we 
compute 1-D spectra of Ee^, Eb^^, and Eb^ vs. k±pi 
by summing the energy of all {kx,ky) modes within a 
given k± shell. Shells are linearly spaced and divided 
into 384 bins; a short-time average over about O.Olr^ 
is performed, as well as an average in z direction. The 
results are displayed in Fig.[T] Here, the solid vertical line 
denotes the boundary to the ’corner modes’, for which the 
angle integration in {kxky) ceases to pick up complete 
circles, causing the artificial spectral break. 

In the range k^Pi ^ 1, an MHD-type spectrum can be 
observed, which exhibits a very small amount of com¬ 
pressive fluctuation energy with a flat spectrum, and 
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Perpendicular wavenumber Pj 


Figure 1. Normalized field energy spectra. Power law expo¬ 
nents obtained from the B± energy spectra within the dotted 
sections are printed into the plot. 


electric and magnetic field energy spectra decaying ap¬ 
proximately with the same power law. The power law 
exponent is close to the Goldreich-Sridhar estimate of - 
5/3 Pl], but the confidence level at small wavenumbers 
is low as there are few modes per shell. 

As the range of k±pi ~ 1 is crossed, all spectra 
steepen, and the turbulence becomes more compressible 
(evidenced by the increased ratio |i?||| / \B±\ ). For 
2 ^ k±pi < 15, all quantities exhibit rather well-defined 
power law spectra, until a further steepening of the spec¬ 
tra sets in at kj_pi « 15, accompanied by a crossing of 
the parallel and perpendicular magnetic fluctuation en¬ 
ergy. These spectral features are consistent with previous 
simulations using a fraction of the present dynamic range 
m- As the choice of parameters is (except for the colli- 
sionality) similar to near-Earth solar wind measurements, 
in Fig. 1 we plot the power law exponent Eg oc ® 
obtained from the measurements of Refs. mm for com¬ 
parison, which agrees within about 15% with our average 
exponent of -3.17, measured between 1 < k±pi < 10. 

Nonlinear energy transfer. In order to study the non¬ 
linear energy transfer, it is useful and necessary to reduce 
the data by subdividing the perpendicular wavenumber 
plane into shells (see also Ref. which we define as 

the region 0 < k± < ko for the 0*** shell and < 

k± < k^TP"!^ for the shells numbered 1 < n < — 1, 

where we set k^ — 0.275 and N = 25. Thus, the entire 
fc_L range present in the simulations is covered, with good 
resolution also for k^pi < 1, while at the same time en¬ 
suring that only the lowest shell 0 < fc_L < /cq contains 
the externally driven modes. 

With this setup, we analyze the net nonlinear shell- 
to-shell energy transfer, which is obtained by summing 
over all q wavenumbers in Eq. (H). The resulting matrix 
(including the symmetric terms due to the antenna, and 
normalized for each k^_ scale) is displayed for the elec- 
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Figure 2. Nonlinear shell-to-shell transfer function for elec¬ 
trons, normalized to the maximum absolute value of each 
wavenumber scale. 


tron species in Fig. Numerical inspection shows that 
the antenna source acts almost exclusively on the lowest 
shell, and diminishes very quickly for higher shell num¬ 
bers. Studying the conservative transfer more closely, 
one can observe that in the range k^pi P 3, while local 
energy transfer dominates, there are some nonlocal con¬ 
tributions connecting disparate k± scales. In the range 
k±Pi > 3, on the other hand, the nonlinear transfer is 
quite local {k± ~ p±), i.e. dominated by direct energy 
transfer between neighboring shells. 

Nonloeal mediation. Beyond the net energy transfer, 
we now extend the analysis to differentiate between dif¬ 
ferent mediators, i.e. q wavenumbers. To this end, we 
evaluate the transfer function of Eq. ([^ with triply fil¬ 
tered inputs, i.e., with fields and distributions condensed 
into shells K,P,Q. Even with the limited number of 
wavenumber shells used here, this diagnostic is extremely 
expensive (approximately oc or about 150,000 core¬ 
hours here), and is thus only evaluated instantaneously 
for a single timestep. Its results can be visualized in a 
compact way, e.g., by means of Kraichnan’s locality func¬ 
tions [5^ . The so-called infrared (IR) locality function is 
defined (following the notation of Ref. [S^) as 


U{kp\kc) 


N 


E 

iC=c+l 


N p p N 

EE + E E 


P=1Q=1 P=lQ=p-|-l 


^K,P,Q 


and retains, for a fixed shell kc with a varying ’probe’ 
wavenumber kp, only transfers for which at least one leg 
p or 5 is smaller than kp. Thus, starting with kp = kc 
(retaining all transfers) and then moving the probe kp 
away from kc, the most local transfers are successively 
removed. For an extensive description of this setup, we 
refer the reader to Sec. V of Ref. [53] . 






















4 



hppi 


Figure 3. Infrared locality functions for several shells fee, nor¬ 
malized to the total nonlinear energy transfer through kc, ver¬ 
sus the probe wavenumber kppi. For the curves with kepi ^ 5, 
a change in slope is apparent when the probe kp crosses the 
ion gyroradius scale. 


For several kc shells, we show the corresponding IR lo¬ 
cality functions 11 (fcp|fcc) /H {kc\kc) in Fig.j^ By plotting 
the curves versus the probe wavenumber kp instead of the 
conventional ratio kp/kc, Fig. [^highlights the existence 
of a meaningful physical scale length at kj_pi ~ 1, indi¬ 
cating a lack of self-similarity. Indeed, the locality func¬ 
tion curves for kePi ^ 5 exhibit a transition in their slope 
that occurs close to the ion gyroradius scale, kpPi ^ 1: 
for kpPi > 1 the nonlinear energy transfer is rather non¬ 
local, with a locality exponent between 2/3 and 1/3; for 
kpPi < 1, a more local exponent of 4/3, as in Navier- 
Stokes turbulence 1221 , is found. As a consequence of 
this property, for 5 < kc ^ 51.2, nonlocal transfers me¬ 
diated by fluctuations in the tail of the MHD range at 
kpPi ^ 1 ^re responsible for at least 30% of the total en¬ 
ergy transfer through these shells. Note that this does 
not contradict the above observation that the net non¬ 
linear transfer for large A:_l is local. Indeed, the nonlinear 
triad fe -t-p -b g = 0 for such nonlocal interactions is char¬ 
acterized by |g| ^ |fe| , IpI and thus |fe| « |p|, consistent 
with a local net transfer between k and p. Finally, we 
note that while all of the above statements were illus¬ 
trated with results for the electron species, the nonlinear 
ion energy transfer (not shown) exhibits the same char¬ 
acteristics, though with an even more pronounced nonlo¬ 
cality (exponent ^ 1/12), and at least 50% of the transfer 
mediated by modes in the tail of the MHD range. 

Collisional dissipation. Next, we study the spectral 
properties of the collisional dissipation rate by measuring 
the contribution of the collision term to the free energy 
balance. The resulting graphs are presented in Fig. 
for both electron and ion species, as well as their sum. 
About 70% of the total dissipation is found to arise from 
electron collisions, which exhibit a broad peak around 
~ 1 ~ 5. Qualitatively, this peak is consistent with 
electron Landau damping acting on the magnetic energy 



Figure 4. Normalized, short-time averaged collisional dissi¬ 
pation for electrons, ions, and its total value. Curves are 
multiplied by k± so the area under the curve is proportional 
to the energy dissipation rate. 

spectrum shown in Figure Despite peaking at these 
relatively small k± wavenumbers, electron dissipation re¬ 
mains strong throughout the spectrum, and begins to in¬ 
tensify somewhat at kj_pi > 30. At k±pi ~ 1, where ion 
transit-time damping is expected to transfer field energy 
to ion particle energy, there is in fact little ion heating. At 
these scales the ion free energy (not shown) is comparable 
to the magnetic fluctuation energy, but it is cascaded to 
smaller scales in both position and velocity space, and is 
dissipated close to the electron gyroradius scale (around 
k±Pi ~ 25). This observation is consistent with an ion 
entropy cascade and the fact that <C i^c niniiiH]. 
Taking into account both species’ contributions, we find 
an essentially flat dissipation spectrum throughout the 
kinetic wavenumber range, contrasting with some inter¬ 
pretations of solar wind data |1[5] which suggested that 
the electron gyroradius scale acts as the dominant dissi¬ 
pation scale. 

Conclusions. In the present study, the first gyrokinetic 
simulation of kinetic Alfven wave turbulence coupling all 
scales from the tail of the MHD range to the electron 
gyroradius scale was performed, with the goal of analyz¬ 
ing fundamental properties of nonlinear energy transfer 
and collisional dissipation for parameters relevant to the 
solar wind. It was found that nonlinear energy transfer 
in the kinetic range, particularly for kj_pi > 5, is con¬ 
siderably more nonlocal than hydrodynamic turbulence, 
as suggested by previous theoretical considerations m, 
and is to a significant percentage (>30%) mediated by 
the tail of the MHD cascade just below k_LPi ~ 1, while 
the net energy transfer occurs mainly between nearest- 
neighbor shells. For Tc/Ti = 1 and /3i = 1, similar to 
the near-Earth solar wind, 70% of the injected energy 
is dissipated through the electron species, whose dissi¬ 
pation spectrum peaks around k±pi ~ 1 — 5, consistent 
with electron Landau damping. The ion free energy, on 
the other hand, is cascaded to small scales and dissipated 
around kj_pi ~ 25. These findings underscore the pres- 
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ence of strong dissipation throughout the kinetic range 
k^_Pi ^ 1, justifying the common notion of a ’dissipa¬ 
tion range’, and demonstrating a coupling across multiple 
scales of both transfer and dissipation. 
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